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Abstract 

Two methods are examined for limit and 
integrator wind-up protection for multi-input, multi- 
output linear controllers subject to actuator constraints. 
The methods begin with an existing linear controller that 
satisfies the specifications for the nominal, small 
perturbation, linear model of the plant. The controllers 
are formulated to include an additional contribution to 
the state derivative calculations. The first method to be 
examined is the multi-variable version of the single- 
input. single-output, high gain. Conventional Anti- 
Windup (CAW) scheme. Except for the actuator limits, 
the CAW scheme is linear. The second scheme to be 
examined, denoted the Modified Anti-Windup (MAW) 
scheme, uses a scalar to modify the magnitude of the 
controller output vector while maintaining the vector 
direction. The calculation of the scalar modifier is a 
nonlinear function of the controller outputs and the 
actuator limits. In both cases the constrained actuator 
is tracked. These two integrator windup protection 
methods are demonstrated on a turbofan engine control 
system with five measurements, four control variables, 
and four actuators. The closed-loop responses of the 
two schemes are compared and contrasted during limit 
operation. The issue of maintaining the direction of the 
controller output vector using the Modified Anti-Windup 
scheme is discussed and the advantages and 
disadvantages of both of the IWP methods are 
presented. 

Nomenclature 

Perturbed vectors 

u c L controller output vector, subject to constraints 
u c controller output vector, (nu x 1) 
u u desired controller output, (nu x 1) 
u, actuator constraints, (nu x I) 
x t . controller state vector, (nx<. x 1) 
x state vector, (nx x 1 ) 
y vector of measured variables, (ny x 1 ) 

z vector of controlled variable, (nzx 1) 


Matrices 

A system matrix of appropriate size 
B input matrix of appropriate size 
C output matrix of appropriate size 
D feedthrough matrix of appropriate size 
L diagonal limit indicator matrix, (nu x nu) 

A CAW Scheme gain matrix, (nx c x nu) 

Scalars 

a variable multiplier for the MAW scheme 

P constant feedback gain for the MAW scheme 

r| constant feedback gain for the CAW scheme 

WF36 fuel flow, lbm/hr 

A8 aft nozzle area, square inches 

Eta ejector butterfly valve angle, degrees 

A78 ventral nozzle area, square inches 

FG9 aft nozzle thrust. Ibf 

FGE ejector thrust, Ibf 

FGV ventral nozzle thrust, Ibf 

N2 fan compressor rotor speed, rpm 

N25 core compressor rotor speed, rpm 

Subscripts 

c controller variable 

d desired value 

i i’th scalar variable in a vector list 

i limit value 

u corresponding to the input vector 

x corresponding to the state vector 

Introduction 

Research in integrated flight and propulsion 
control performed at the Advanced Control Technology 
Branch of the NASA Lewis Research Center has 
resulted in an Integrated Method for Propulsion and 
Airframe Control (IMPAC) design [ I]. To evaluate this 
method, an example control system is being designed 
for a model of a Short Take-Off and Vertical Landing 
(STOVL) aircraft during transition flight. The steps of 
the IMPAC design are shown in Figure I . Step five of 
the method is the nonlinear control design for the 
various subsystems. The nonlinear control design 
includes the addition of washout filters and trim logic, 
gain scheduling, limit logic, and limit and integrator 
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windup protection (IWP). Preliminary evaluations of 
the NASA Lewis linear integrated control system for 
large command inputs have indicated the need for limit 
and integrator windup protection on the propulsion and 
aircraft subsystems. Integrator windup protection is 
common in control implementations, and schemes for 
single-input, single-output (S1SO) systems have been 
discussed in the literature. General, multi-input, multi- 
output (MIMO) integrator wind-up protection methods 
have not received a great deal of attention in the 
literature, until recently. This paper examines and 
compares two of the available multi-variable integrator 
windup protection schemes for linear control systems 
with constraints imposed upon the range of the 
actuators. 
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Figure 1 IMP AC Methodology Flowchart 


In the following, the integrator wind-up and 
limit protection problem for general, multi-input, multi- 
output, linear controllers is presented. Then, the 
literature on multi-variable integrator windup methods 
is briefly reviewed. Next, two methods, the 
Conventional Anti-Windup (CAW) and the Modified 
Anti-Windup (MAW) schemes are examined in detail. 
To test these two schemes, an example problem is 
described, consisting of a linear control system for a 
linear model of a turbofan engine with actuator range 
limits. The closed-loop performance of the CAW and 
the MAW schemes during limit operation are compared. 
Finally, the advantages and disadvantages of these two 
IWP methods are discussed. 



Figure 2 Typical Limit Protection Problem Structure 


The Integrator Windup Problem 

All control implementations have the possibility 
of encountering constraints such as actuator range and 
rate limits or system safety limits reflected back upon 
the plant inputs. When limits occur, the control system 
integrals can wind up if the controller was designed to 
drive the steady-state errors to zero. The windup 
appears as an increase in the magnitude of the 
integrators beyond that which would normally be 
observed during linear operation for the same inputs. A 
simple example of integrator windup is an error driven 
SISO proportional plus integral plus derivative, (PID), 
control system with a limited actuator. With the 
actuator constrained, the error can not be driven to zero. 
As long as the error is nonzero and of constant 
numerical sign, the integral term in the controller will 
increase in magnitude. Integrator windup is only 
unbounded for a constant constraint. Transient 
encounters with the actuator constraints cause short 
duration integrator windup that can lead to oscillations 
such as limit cycles or delayed responses due to the 
controller holding the limit value until the integrator 
"unwinds". A typical structure for the limit protection 
problem is shown in Figure 2. The structure consists of 
a plant, G(s), a nominal linear controller, K(s), and 
actuator range limits. The linear, time invariant, small 
perturbation, model of the plant, G(s), is represented in 
state space format as follows: 

x = Ax + Bu 


z 

k 


= Cx + Du 


( 2 ) 


The linear control system, K(s), satisfies the nominal 
closed-loop, small perturbation specifications and is 
represented by the following state-space structure: 
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(3) 


Integrator Windup Protection Schemes 


= A c x e + B c 
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C x 


e 

ly. 


(4) 


Figure 2 shows that the controller output vector, is 
the input to the limit block. The vector output of the 
limit block is u c L . As long as the magnitude of the Tth 
element of u c , u d , is less than the i’th actuator limit 
value, u^, then u c L , = u d , otherwise u d L is equal to the 
corresponding actuator limit value, The actuator 
limit error, defined as e n =u c L -u c , is fed back to the IWP 
scheme. Actuator limits are typically taken into 
consideration during the selection of the input scale 
factors and weights that are used in the nominal control 
design. The actuator limits will be encountered 
whenever maximum performance is required as the 
constraints define the maximum system performance. 
IWP will be required during this limit operation. The 
specifications on rWP schemes, as discussed in 
Reference [2], are as follows: 

1) An actuator limit represents a constraint in the input 
space and must be observable. 

2) The IWP scheme should not contribute to the 
nominal control system when the system is 
unconstrained. This indicates that the IWP scheme 
should be memoryless. 

3) The controller with IWP should track the limited 
actuator value such that the transition from limit 
mode to unlimited mode is smooth, (a "bumpless 
transfer"). 

4) The IWP scheme should be closed-loop stable for 
all possible limit combinations. 

5) When actuator limits occur, the IWP scheme should 
maintain system performance. If it is not possible 
to maintain the closed-loop performance, then the 
IWP scheme should degrade the system 
performance gracefully, while maintaining stability. 

In the following, the actuator limit values are assumed 
to be known so that the controller outputs can be limited 
within the control itself, before the actuator command 
goes to the actuator. Also, the IWP schemes discussed 
here are memoryless. Guarantees for closed-loop 
system stability with actuator constraints and guarantees 
for closed-loop performance for systems with redundant 
actuators are both current research topic and are outside 
of the scope of this paper. The literature on multi-input, 
multi-output limit and integrator windup protection 
schemes will now be reviewed. 


The literature on limit and integrator windup 
protection discusses several techniques. For SISO 
systems, the error that drives the integral error term is 
typically zeroed when the system encounters a limit. 
Other SISO schemes [3], have fed back the observed 
actuator error, e B , as done in Figure 2. This error is 
used to drive the controller to track the limited actuator, 
putting a limit tracking proportional regulator upon the 
controller itself. Several versions of this e Q feedback 
approach have been extended to Ml MO systems and 
they are known, collectively, as the Conventional Anti- 
Windup (CAW) method. In a M1MO system, an 
encounter with a single limit changes the direction of 
controller output vector, u c . Figure 3 shows such an 
example for a two dimensional vector during an 
encounter with the limit for the first actuator such that 
u dl =UL,. The scalar approach modifies the magnitude 
and not the direction of u c , where as the CAW scheme 
truncates the limited actuator value. Some systems are 
sensitive to changes in the direction of u c and for these 
systems an actuator limit can result in poor closed-loop 
performance. 



Figure 3 2D Actuator Directionality 


Reference [4] addresses this directionality issue for 
stable plants and square controllers by multiplying the 
error vector, e, by a scalar, A(t), (0 < A £ 1). This 
"error governor" is a function of the boundary of the 
admissible controller state which can be calculated by 
reflecting the actuator limits onto the controller states 
and by predicting bounds for the controller state 
trajectory. The scalar multiplier prevents the controller 
from exceeding the actuator limit and zeros the error 
input to the controller when a saturation is reached. 
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The usefulness of this method is limited because it 
requires the solution to C c e A '\-u, for the controller 
state, x a such that any possible combination of the 
actuator limits, u,, are not violated in order to determine 
the region for the admissible controller state. This 
problem is formidable, particularly if one considers a 
controller with time varying limits, or a gain scheduled 
controller, or nonsymmetric bounds for u,. Also, the 
solution to this problem is not unique if there are more 
controller state variables than controller outputs. 

Reference [5] presents an alternative approach to 
the limit problem using linear robust control design 
tools. In Reference [5], the saturation function is re- 
formulated as a dead-zone nonlinearity to model the 
actuator limits as a sector bounded uncertainty. Then, 
linear multi-variable control design methods guarantee 
the system closed-loop robust stability for specific levels 
of this "uncertainty". Adding actuator limit uncertainty 
to all of the plant inputs yields a design that is similar 
to one obtained by increasing the weights on all inputs 
and does not directly provide IWP. CAW limit 
protection can be built into this design plant structure if 
a zero steady-state error specification is provided on the 
actuator error, e^, and if this error is observed by the 
controller. This approach combines the design of the 
nominal, linear controller and the IWP into one step and 
attempts to simultaneously satisfy the specification of 
both the nominal controller and the IWP. This approach 
results in a trade-off of the small perturbation, nominal, 
closed-loop performance for stability robustness to limit 
operation and may be conservative. 

Reference [2] presents a general theory of anti- 
windup/bumpless transfer methods and uses this theory 
to compare the structure of several versions of the 
Conventional Anti-Windup (CAW) method. A version 
of the CAW method is demonstrated in Reference [6], 
Reference [2] alludes to an IWP method denoted as the 
Modified Anti-Windup (MAW) method and this method 
was demonstrated in Reference [7] on a academic, two 
input, two state, two output controller. The MAW 
scheme uses a scalar multiplier similar to that in 
Reference [4], except now the scalar multiplier is 
applied to the controller output instead of the input. 
Both the CAW and the MAW schemes will be reviewed 
in the following. 


The Conventional Anti-Windup Scheme 

There are several variations on the CAW scheme 
as discussed in reference [2]. Figure 4 shows the 
structure for the version of the CAW scheme that is 
implemented in this paper, which includes the matrices 
for the nominal controller, [A* C c D c ], from 


H 



Figure 4 Implementation of the CAW Scheme 


equations (3) and (4), plus an additional contribution to 
the controller state derivative calculation, Ae 0 , as shown 
in Figure 4. This term is memoryless, since A is a 
constant gain matrix. Using e B =u c L -u c , u f as defined in 
equation (4), and the fact that u rt L =u u whenever the fth 
actuator is limited, results in the following CAW 
controller state equation: 


*c = Vc + B « 


ALe„ 


= (A e -ALC c )x e +(B c - ALDj 


(5) 


+ALu, 


Note in equation (5), u, is a vector of actuator limits. 
The CAW controller output equation is the same as 
equation (4). Note that the nominal controller is active 
as long as L is zero. L is a diagonal limit indicator 
matrix of Ts and 0’s, one indicating a limited actuator. 
For the case where the third actuator is limited for a 
controller with 4 outputs, L is defined below as: 


L « 


0 0 0 0 
0 0 0 0 
0 0 10 
0 0 0 0 


( 6 ) 


The controller stability can be examined by studying the 
finite number of combinations of actuator limits using 
(A c -ALC e ). The matrix A is the only design parameter 
for the CAW scheme. The approach taken here for 
designing A follows the "tracking mode" implementation 
suggested by AstrOm and Wittenmark, as discussed in 
Reference [2]. For this implementation the nominal 
controller is in modal form and A was designed to 
contribute to the diagonal of controller system matrix, 
A^ This is accomplished if ALC is diagonal. This is 
approximated using A m r\C c ( # indicates the 
pseudoinverse and n is a scalar). All possible limit 
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combination were then checked for controller stability. 
There are other methods for selecting A and this 
explains the family of CAW schemes. The idea behind 
this CAW scheme is to track any actuator limit by back- 
calculating the controller state such that the actuator 
limits are satisfied. 


= (A c + P(«-l)0x c ♦ B c 


|e 

ly. 


( 8 ) 


= C c\ 


aD„ 


( 9 ) 
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Figure 5 Implementation of the MAW Scheme 
The Modified Anti-Windup Scheme 


The MAW scheme provides IWP without 
modifying the direction of the control calculated output 
vector, as shown in Figure 3. Figure 5 shows the 
implementation structure of the MAW scheme that is 
used in this paper. Note that the nominal controller, as 
described in equations (3) and (4), is active as long as 
a=l. When a*l, the nominal controller is modified by 
the two additional blocks shown in Figure 5. One of 
these blocks scales the feedthrough term D c [e T y T ] T by 
the scalar, a. The other block "effectively" scales the 
C c \ term by a using additional state feedback to the 
controller state derivative calculation using a gain of 
p(a-l). P is a constant scalar design parameter and is 
always greater than 0. a is defined as follows: 


a(t) = 1, if u r \ = u c , , else 


a(t) « min 
i 


u c(t), 


U e (t), * 0 


(7) 


The above ratio is the result of the division of the fth 
element of output from the limit block, u ri L , by the i’th 
element of controller output, u d . Note that this is the 
ratio of perturbed variables and it has been assumed that 
u c and u c L are of the same numerical sign. Whenever 
|u d |>|u u |, for any i, then a is equal to the smallest 
of the ratios indicated in equation 7. This implies 
0<a<I. The controller equations for the MAW scheme 
corresponding to Figure 5 are as follows: 


The MAW structure only modifies the real part of the 
controller eigenvalues, regardless of the controller 
coordinates. This can be seen by applying a coordinate 
transformation matrix, T, to the MAW controller system 
matrix, (A c +p(a-l)I), as shown below: 

T-'(A+p(a-l)l)T = T "'AT + p(a-l)l (10) 

T does not change the scalar term and thus the MAW 
scheme is independent of the controller coordinates. 
Also, since p>0 and (ot-l)<0 for 0<a<l, the MAW 
scheme moves the real part of the eigenvalues of the 
nominal controller further into the left half plane by the 
amount P(a-l). Thus the MAW can never make the 
controller unstable. It would be ideal if this method 
guaranteed the stability of the closed-loop system, but 
this is an area for future study. 


scalar 



Figure 6 Idea Behind the MAW Scheme 

The idea behind the MAW scheme is shown in 
Figure 6 by replacing the limit block with the scalar 
variable, a. Note that a scales all of the actuator values 
when any actuator is constrained. The desired control 
calculated output value is defined as u u =au r where a is 
defined in equation (7). This scalar multiplier block 
does not provide IWP, so a is moved into the control 
system block using block diagram manipulation as 
shown in Figure 7. The feedthrough term, D c [e J y 1 ] 1 
can be scaled directly by a, but the controller state must 
by modified through the controller state derivative, 
unless a resctable integrator is used. There are many 
ways to modify the controller state derivative and the 
following scheme was selected because it provides fast 
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limit tracking and because it resembles the CAW 
scheme. By assigning a portion of the actuator error, e*, 
to the controller state and defining the desired state 
vector, the state error vector can now be 

defined as e^x^-x^ax^-x^a-l )x*. The desired output 
equation can now be written as follows: 


= C c ox e * aD c 

" C e X cd + aD c 


(ID 


A proportional gain controller acting on the state error 
vector can be used to drive x c to x^. Figure 5 shows 
the resulting structure with p as the proportional gain. 
It is interesting to note that the MAW scheme is 
independent of the controller matrices and thus would 
not have to be scheduled if the controller were gain 
scheduled. 



Figure 7 Obtaining IWP for the MAW Scheme 
via Block Diagram Manipulation 


Turbofan Engine Limit Model 

Reference [6] discusses the implementation and 
response of a closed-loop engine control system as part 
of the evaluation of an integrated flight/propulsion 
control system. A linear version of the closed-loop 
propulsion control system and plant used in Reference 
[6] is used in the following example. Figure 8 shows 
the closed-loop system, comprised of a nominal 
controller with IWP, linear actuator models, actuator 
limits, and a linear engine model in the form given by 



Figure 8 Closed-loop Turbo fan Engine Example 

equations (1) and (2). The inputs, controlled variables, 
extra feedback, and state are defined as follows: 

u = [WF36, A8, Eta. A78] T 

z = [N2, FG9, FGE, FGV] T (12) 

y = N2 

x = [N2, N25 ] t 

The above variables are defined in the nomenclature list. 
The limits for this example are hard actuator limits 
reflected onto the linear perturbation variables and these 
limits were further reduced to study the limit operations 
of the IWP schemes. For example, at the design point 
for this linear engine model, the nominal value for the 
ejector butterfly valve angle, Eta, was 64 degrees and 
the maximum butterfly valve angle is 90 degrees (full 
open). This yields a perturbation limit of 26 degrees for 
Eta, but the perturbation limit was artificially lowered to 
9 degrees to study the limit operation of the two IWP 
schemes. Other perturbation limits were similarly 
modified during this study. The nominal control system 
design is a linear, perturbation design as discussed in 
Reference [8]. The engine controller is 7th order. The 
reason for this relatively high order controller, (the plant 
is second order), is that the original controller was 
designed with a higher order model of the engine which 
included temperature and pressure state variables. The 
eigenvalues for the open-loop engine, the nominal 
controller and the controller with MAW active are 
shown in Table 1. The MAW controller eigenvalues 
were calculated using f5-10 and a=0.9 for this example 
but a variations in the range from 0.8 to 1.0 were 
observed in the simulations. The MAW scheme makes 
the real part of all the eigenvalues more negative by the 
amount p(a-l), as seen in Table I. Table 2 shows the 
controller eigenvalues with the CAW scheme for various 
combinations of actuator limits. For example, the 
matrix L-diag[0 0 1 0] corresponds to the case with the 
ejector butterfly valve. Eta, is limited, n* 10 is used in 
this example to provide reasonable limit tracking. The 
closed-loop results for the CAW and MAW controller 
schemes will now be examined. 
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Results of the Turbofan Engine Example 

Various commands were simulated for the closed- 
loop system using both IWP schemes and the results 
were compared. One response that gives a good 
summary of the results will be examined. The 
command is an ejector thrust command, FGE^. The 
command of FGE alone is atypical, but it allows the 
changes in coupling to be more easily seen since the 
other loops are trying to regulate to zero. This single 
command does not detract from the comparison of the 
IWP schemes. The FGE^ transiently pushes the Eta 
actuator into the limit. Figure 9 shows the controller 
outputs, u, for the three cases: (1) the Modified 
AntiWindup (MAW) scheme; (2) the Conventional Anti- 
Windup scheme; (3) the unlimited case. Note that ETA 


equal to nine is the only limit encountered. The last 
plot of Figure 9, the controller output value, ETA_c, 
shows that both IWP schemes closely tracks the limit 
value. The response for WF36, A8, and A78 in the 
CAW case are similar to the response for the unlimited 
case, but for the MAW scheme, all the actuators are all 
scaled back by a. Figure 10 compares the closed-loop 
system outputs. FGE command is also shown in the 
third plot of Figure 10 for reference. Note that the level 
of coupling for the CAW scheme appears to be slightly 
greater that the coupling for the case with the MAW 
scheme. However, the FGE command tracking for the 
controller with CAW appears to be better than the 
command tracking for the controller with the MAW 
scheme, oc is shown in Figure 10 for reference. 
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Figure 9 Comparison of Actuator Responses Figure 10 Comparison of Closed-Loop Outputs 
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From Figure 9 and Figure 10 it is difficult to determine 
if one method is superior to the other and that was 
fairly typical of all the responses that were examined. 
This turbofan engine example does not reveal a clearly 
superior method, perhaps because this example is not 
particularly sensitive to changes in the direction of the 
control vector. Another example that does exhibit 
sensitivity to the direction of the control vector is the 
aircraft pitch control problem. In the STOVL aircraft 
example currently being studied [6], FGE and FGV 
from the current turbofan engine example balance the 
pitching moment to provide pitch control for a STOVL 
aircraft. Preliminary results from this study indicate a 
distinct advantage for the MAW over the CAW scheme 
and these results will be published in the future. 

One claim of the scalar, MAW scheme is that it 
does not modify the direction of the input vector. This 
needs clarification. Multiplying a vector by oc(l^) does 
not change the vector direction at that time, t„. This 
does not mean the vector direction of the limited system 
with MAW will be a scalar multiple of the unlimited 
closed-loop system. For the vector direction, relative to 
the unlimited case, to be maintained exactly, the zero 
crossings for the unlimited and the MAW schemes 
would have to be identical. To probe this directionality 
question further, consider the constant value of a shown 
in Figure 1 la. In this case, the constant cc scales the 
actuator and the inverse of <x scales the corresponding 
measurements. This has the effect of fictitiously 
making the plant look larger than it is to the controller, 
(y/a, 0<a<l). Time history responses of this closed 
loop system for various constant values of a are all 
properly scaled and the zero crossings of the controller 
outputs and plant outputs are identical to the case for 
a=l. In fact, block diagram manipulation shows that 
for a=constant, this is really just a scalar reduction of 
the command as shown in Figure I lb, ( 1/a K(s)a = 
K(s) ). Consider a time varying value for a, a(t), in 
Figure 11a. For this case, directionality relative to the 
case where a=constant, is not maintained. The reason 
for this is that the plant has memory and the plant state 
vector does not get appropriately scaled by a time 
varying scalar. But Figure lib does help to understand 
how the MAW works. The MAW scheme will appear 
as a variable rate limit on all the commands. There still 
remains a question regarding the affect the MAW has 
on the closed-loop feedback properties. 

Several practical aspects of CAW and the MAW 
schemes need to be addressed. One of the advantages 
of the MAW scheme is that it does not have to be 
altered if the controller is gain scheduled. This is an 
advantage over the CAW scheme where the 


pseudoinverse of the C r matrix would have to be 
calculated online as it was scheduled throughout the 
operating range. One disadvantage of the MAW 
scheme is that it is implemented on the perturbation 
variables. This requires that the numeric signs of u c and 
u r L must be the same, which may not always be true 
depending on the cause of the actuator limit. For 
example engine safety limits like the fan surge margin 
may have internal logic that makes it possible to have 
the perturbed variables u c L and u c of different numeric 
signs. Additional logic would be required to handle 
these special cases. The CAW scheme is based on the 
actuator error, e g , and does not have this problem. In 
terms of stability, it is desired to have a test for closed- 
loop stability. While this is not yet available, the MAW 
scheme can not make the controller unstable, but on the 
other hand the CAW scheme has a finite number of 
cases to test for closed-loop stability. 



Figure 1 la Scaling the Plant 


Kfs) 



Figure lib Scalar Modification of Commands 
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Conclusions 

The conventional anti-windup (CAW) and the 
modified anti-wlndup (MAW), multi-variable limit and 
integrator windup protection (IWP) schemes were 
described and compared. A design example consisting 
of a model of a turbofan engine control system was 
presented to compare the closed-loop limit operation of 
the two IWP schemes. The example system presented 
does not exhibit sensitivity to changes in the direction 
of the control vector and thus neither of the IWP 
schemes were shown to be clearly superior. Both 
schemes were successful in providing integrator windup 
protection and limit tracking. Nevertheless, the MAW 
does look promising as it attempts to maintain the loop 
properties by moving the limit problem out of the 
feedback loop and on to the commands as rate limits. 
An example system that does exhibit sensitivity to 
changes in the direction of the control vector, the pitch 
control of a STOVL aircraft, is currently being 
investigated. 
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Table 1 Eigenvalues for the 
Open- loop Nominal Controller 

Plant Controller with MAW 

-2.97 

-0.004 

-1.004 

-6,33 

-0.007 

-1.007 


-0.011 

-1.011 


-0.011 

-1.011 


-0.011 

-1.011 


-0.171 

-1.171 


-2.896 

-3 . 896 



Table 

2 


Eigenvalues for CAW Scheme 

L- [0100] 

[0111] 

[1010] 

[1111] 

-0.004 

-0.004 

-0.004 

-0.004 

-0.008 

-0.009 

-0.010 

-0.010 

-0.011 

-0.150 

-0.011 

-2.285 

-0.011 

-2.260 

-0.168 

-20.01 

-0.011 

-20.01 

-2.781 

-20.01 

-2.801 

-20.01 

-20.01 

-20.01 

-20.11 

-20.64 

-20.12 

-20.64 


L - [0100] Is a diagonal matrix 

with 0, 1, 0, 0 on the diagonal 
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